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Abstract 

For energies relevant to future linear colliders, y/s > 500 GeV, the WW fusion 
channel dominates the Higgs boson production cross section e + e~ — > uvh°. We have 
calculated the one-loop corrections to this process due to fermion and sfermion loops 
in the context of the MSSM. As a special case, the contribution of the fermion loops 
in the SM has also been studied. In general, the correction is negative and sizeable 
of the order of 10%, the bulk of it being due to fermion loops. 



As no Higgs boson could be detected so far, the search for the Higgs boson as the 
primary goal of high energy physics will continue. The four LEP experiments delivered 
a lower bound for the Standard Model (SM) Higgs mass, m h > 114 GeV 0. In e + e~ 
collisions, for energies > 200 GeV, the production of a single Higgs boson plus missing 
energy starts to be dominated by WW fusion [3|, ||, that is e + e~ — > v e v e WW — > v e v e h°, 
whereas the Higgsstrahlung process || e + e~ — > Zh° — > vv h becomes less important, h 
denotes the SM Higgs boson or the light Higgs boson of the Minimal Supersymmetric 
Standard Model (MSSM), while H° is the heavy CP-even Higgs boson in the MSSM. 
The rates for the ZZ fusion are generally one order of magnitude smaller than those of 
the WW channel. 

At LHC, in pp collisions, the gluon-gluon fusion mechanism provides the dominant 
contribution to Higgs production. The next important Higgs production channel is vector- 
boson fusion VV —>■ h°/H°. In particular, it provides an additional event signature 
due to the two energetic forward jets. Recently, it has been argued that the channels 
WW — > h°/H° — > rf and WW, can serve as suitable search channels at LHC even for a 
Higgs boson mass of mh ~ 120 GeV ||. 

At Tevatron, with pp collisions at 2 TeV, the WW fusion process plays a less im- 
portant role. The g g fusion is the dominant process for Higgs production there, but the 
Higgsstrahlung q q —>■ W —>■ Wh° is larger than the WW fusion for < 180 GeV 0. 

It is worth mentioning that WW (and ZZ) fusion is also the most important Higgs 
boson production mechanism in ep collisions. Actually, it was shown || that the ep- 
option at LHC would offer the best opportunity to search for a Higgs boson in the mass 
range < 140 GeV. 

In this paper, we have calculated the leading one-loop corrections to the WWh° vertex 
in the MSSM by taking into account fermion/sfermion loops. They are supposed to be 
the dominant corrections due to the Yukawa couplings. We have applied them to the 
Higgs boson production in e + e~ annihilation in the energy range \/s = 0.5 — 3 TeV, i. e. 
e + e~ — > v e v e WW — > z/ e i/ e /i°. We included the Higgsstrahlung process e + e _ — > Zh° — > 
vv h° and the interference between these two mechanisms. Because the Higgsstrahlung 
process is much smaller in this range, we have neglected its radiative corrections. We have 
also discussed the SM case where we have studied the dependence on the Higgs boson 
mass. As to the one-loop corrections to the WWh° vertex in the SM, we quite generally 
refer to the review article 11. For this coupling also QCD corrections were included, the 
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As for energies a/s > 500 GeV the dominant channel e + e — > vvh° is by far the WW 
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Figure 1: The Feynman graphs for the process e + e — > uuh°. Note that for lA^ 66 ] 2 one 
has to sum over all three neutrino flavors. 

fusion, in the context of our calculation only corrections to the WWh° vertex must be 
evaluated, see Fig. At the one-loop level the Lagrangian for the WWh° coupling can 
be written as 

£=(f&w9r+(*$wY) h°W+W~, (1) 

where the tree-level coupling in the MSSM can be read as g^w = 9 m w sin(/3 — a). The 
one-loop part of the WWh° coupling in Eq. (JJ) can be expressed in terms of all possible 
form factors^] as 

(A^)"" = F 00 g^ + F n k$kl + F 22 k%k» + F 12 k<tk u 2 + F 2X k^k\ + 1 F e e^ s k lp k 2S , (2) 

where k\^ denote the four-momenta of the off-shell P^-bosons. At tree-level only the 
structure with g^ v is present, and therefore all form factors but F 00 have to be ultra 
violet (UV) finite without being renormalized. By adding appropriate counter terms also 
the form factor F 00 is rendered UV finite. For the renormalization procedure the on-shell 
scheme has been adopted. 

Actually, for the calculation of the one-loop corrected WWh° vertex we compute the 
vertex and the wave-function corrections stemming from the graphs of Fig. ||], as well as 
the coupling correction, 

( a h° fx h a (v)\fJ.is . (, h°(w) . r h°(c)\ uv (n\ 

[Ag ww J ={Sg W w ) + \ 6 9ww + s 9ww ■ ) 9 ■ ( 3 ) 

Comparing Eq. (Q) with Eq. ([|), it is obvious that the counter terms, i. e. the wave- 
function and coupling corrections, contribute only to F 00 . The wave-function correction 
is 

^9ww^ = 9ww [\ {$ Z H)h°h° + SZ W ) + | g^ w (8Z H ) H o h o , (4) 
1 As the expressions for F 00 , F n ,...,F e are quite lengthy they will be presented elsewhere fl2|| . 
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Figure 2: The Feynman graphs that contribute to the vertex corrections defined in Eq. 
and the wave-function in Eq. (|J). / (/') denotes the up (down) type fermion. 



where 8Zh and 5Zw are the symmetrized Higgs boson and the VT-boson wave-function 
correction, respectively, calculated from the corresponding graphs of Fig. In the case 
of the off-shell P^-bosons, 5Zw has the form 
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The coupling correction is 

5a h ° (c) - 
°9ww ~ 



5g 5mw\ h° sin 2/3 5tan/3 H o 

I 1 9ww "I 7, 7 7T 9ww 

g m w J 2 tanp 



(5) 



(6) 



The expressions on the right-hand sides of Eqs. (§-§D can be found in Ref . [|13| . Especially, 
we fixed the counter term 5tan/5 by the on-shell condition Imll^z(m^) = 0, where 
^Az{ m2 A> is the renormalized self-energy for the mixing of the pseudo-scalar Higgs boson 
A and Z-boson. 
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Having calculated the form factors of Eq. (fj), one can proceed to the calculation of 
the one-loop corrected cross section. The tree-level part of the amplitude was already 
calculated in Ref. Q. Here we will include also the one-loop correction, that is the 
interferences between the tree-level and the one-loop amplitudes of Fig. For these 
terms we get 
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AM w (M^ cc y] = ±g$ w g 4 (2F 00 Pl -p 4P2 - P3 + F 21 S)l[ 



1 



and 



23ft 



AM W (Mf e y 



where 



=1,2 



m 



(7) 



W) 



9' 



cos 2 6 



w 



(2F 00 Pl -p 4 p 2 -p 3 + F 21 S) 



x 



4 



(s? - m%) {{ql - mf ) 2 + m| V%) II k 



m 



w 



S = {pi ■ Pa+ P2 ■ P3) {Pi ■ P2P3 ■ Pa + Pi ■ P4P2 ■ P3 - Pi ■ P3P2 ■ Pa) 

- 2(p 1 -p 2 +p 3 -p 4 )p 1 -pAP2-P3, 



(9) 



T z is the total Z-boson width, and C£ = — | + sin 2 ^vK- The kinematics has been chosen 
as e~(pi) e + (p 2 ) — > v e {pz) v e {j>A) h°(p). We have also defined k\ = p 3 — pi, k 2 = p 2 — p 4 , 
qi=Pi+ P2 and q 2 =p 3 + p 4 . 

For the calculation of the cross section at tree-level, it is possible to perform some of 
the phase space integrations analytically and the rest of them numerically || f| . However 
including the one-loop correction terms of Eq. (|7|) and (^), it is impossible to perform 
any of these integrations analytically. Therefore, we have performed the integrations 
numerically, using efficient numerical integration subroutines found in the NAG library. 
We have also checked that for the tree-level case our completely numerical calculation 
agrees with high accuracy with the semi-analytical results of Ref. || . 

Before embarking on discussing our numerical findings, we will make some comments 
concerning some details of our calculation. For the calculation of the supersymmetric 
(SUSY) Higgs boson masses and the Higgs mixing angle a, a computer programme based 
on Ref. has been used. The tree-level WWh° coupling for values of tan f3 > 5 as 
preferred by the LEP Higgs boson searches, mimics the SM one. For the calculation of 



2 As it has already been mentioned the corresponding corrections to the Higgsstrahlung process can 
be safely ignored due to absolute dominance of the WW fusion for y/s > 500 GeV. 
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Figure 3: The various cross sections as functions of yfs (left). The clotted-dashed line 
represents the tree-level cross section a^ w , the dotted line a^ w + <jQ~ str . The dashed 
line includes also the interference term o"o nterf ' and represents the total tree-level cross 
section. The solid line includes the one-loop correction. The SUSY parameters are: 
tan/5 = 10, fi = -100 GeV, A = -500 GeV, ttiq = 300 GeV, M A = 500 GeV, and 
M 2 = 400 GeV. For the same set of parameters also the tree-level cross section <jq and 
the one-loop corrected a are plotted for yfs up to 3 TeV (right). 

the fermion/sfermion one-loop corrections to the WWh° vertex, the contribution of the 
third family of fermions/sfermions has been taken into account. This contribution turns 
out to be the dominant one, in comparison with the first two families corrections, due 
to the large values of the Yukawa couplings ht and h^. The effect of the running of the 
coupling constants g and g' has been taken into account. 

We now turn to the discussion of our numerical results. In Fig. |3| (left) we can see 
the various cross sections as a function of yfs for values up to 1 TeV. The dotted-dashed 
line represents the contribution from the WW channel at tree-level alone, whereas the 
dotted line includes the Higgsstrahlung contribution as well. The dashed line comprises 
in addition the interference between the WW channel and Higgsstrahlung. One can 
perceive that the size of this interference term is extremely small, and for this reason 
the difference between the dotted and dashed lines is rather minute. It is also clear that 
for yfs > 500 GeV the WW fusion contribution dominates the total cross section for 
the Higgs production e + e~ — > uvh°. Actually, for yfs > 800 GeV the total tree-level 
cross section is due to WW fusion. In the solid line we have taken into account the 
one- loop correction from the fermion/sfermion loops. The correction is always negative 
with a significant size of the order of 10%. This can also be seen from Fig. |3] (right), 
where we have plotted the tree-level cross section <jq (dashed line) and the one-loop 
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Figure 4: The relative correction Act/cto as a function of \/s (Act = ct — o~o), where cto 
is the tree-level and ct the one-loop corrected cross section. The solid and dashed lines 
correspond to two different choices of the SUSY parameters, as described in the text. 

corrected cross section ct (solid line) for energies up to 3 TeV. For simplicity, for all 
plots we have used A t = At, = A T = A, {m^, m^,m^,m^} = ^, 1, 1} uiq and 
Mi = §M 2 tan 2 # w . The choice of a common trilinear coupling and the correlation 
between the soft sfermion masses are inspired by unification. For the plots in Fig. |3] we 
have taken: tan/3 = 10, fj, = -100 GeV, A = -500 GeV, m Q = 300 GeV, M A = 500 GeV, 
and M.2 = 400 GeV. Choosing different sets of parameters, the basic characteristics of 
these plots remain indifferent. In fact, the soft gaugino masses M 12 affect only the Higgs 
sector through radiative corrections. 

In Fig. f| the relative correction Act/ ct is presented as a function of ^Js for two different 
sets of parameters. The solid line corresponds to the set tan/3 = 10, n — —100 GeV, 
A = -500 GeV, tuq = 300 GeV, M A = 500 GeV, and M 2 = 400 GeV, whereas for the 
dashed line we have taken tan/3 = 40, fi = —300 GeV and A = —100 GeV, keeping the 
rest of them unchanged. This figure shows that the size of the one-loop correction to the 
Higgs production cross section is practically constant for y^i > 500 GeV and weighs about 
—15%, almost independently of the choice of the SUSY parameters. The reason for this 
weak dependence is that the one-loop corrections are dominated by the fermion loops, and 
therefore the total correction is not very sensitive to the choice of the SUSY parameters. 
In this region, a/s > 500 GeV, one can compare this constant correction with the effective 
approximation of Ref. [HJ, where one only corrects the WWh° coupling. Although the 
sign of this approximation is correct, it does not fully account for the whole effect. 

In Fig. ^ we have plotted the cross section as a function of uih for the SM case, 
for y^s = 0.8 TeV (red lines) and 1 TeV (black lines). The dashed lines correspond to 
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Figure 5: Cross sections for the SM case, for \fs = 0.8 TeV (red lines) and 1 TeV (black 
lines). The dashed lines correspond to the tree-level cross section, whereas the solid lines 
to the one-loop corrected one. 

the tree-level cross section for e + e~ —* uuh°, whereas the solid lines contains the one- 
loop correction stemming from the fermion loops. In addition, the couplings have been 
adjusted to the SM corresponding couplings. The plot exhibits the expected dependence 
of the cross section on m^. What must be noticed is that especially for small Higgs boson 
masses < 200 GeV, the size of the fermion loops correction becomes important for the 
correct determination of the Higgs boson mass. 

Finally, Fig. |6| exhibits the percentage of the sfermion loops to the total one-loop 
correction as a function of tan/3 (left) and \i (right), for two different values of /i and 
tan (3, respectively, as shown in the figure. In the left (right) figure we have chosen 
A = -100 GeV (A = -400 GeV). The rest of SUSY parameters are: m Q = 300 GeV, 
Ma = 500 GeV, and M 2 = 400 GeV. Here -^s has been fixed to 1 TeV. The grey area in 
the right figure is excluded due to the chargino mass bound. We see that the maximum 
value of order 10% can be achieved for large values of fi and tan/3. There, due to the 
significant mixing in the stop and sbottom sector, the contribution of stops and sbottoms 
in the loops is enhanced. For these values of SUSY parameters the sfermion masses 
approach their experimental lower bounds. But even there the dominant correction, at 
least 90% of the total correction, is due to the fermion loops. 

In conclusion, we have calculated the fermion/sfermion loops corrections to the single 
Higgs boson production e + e~ — > uuh in the context of the MSSM and SM. They are 
supposed to be the dominant radiative corrections. For energies relevant to the future 
linear colliders, a/s > 500 GeV, the WW fusion channel dominates the cross section. In 
general, the correction due to fermion/sfermion loops is negative and yields a correction 
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Figure 6: The percentage of the sfermions to the total one-loop correction as a function of 
tan/? (left) and fi (right). The rest of the SUSY parameters have been fixed as described 
in the text. Here \/s = 1 TeV. The grey area in the right figure is excluded due to the 
chargino mass bound. 

to the cross section of the order of —10%. The bulk of this correction stems from the 
fermion loops, and usually turns to be more than 90% of the total correction. For the 
case of maximal mixing in the sfermion mass matrices, the contribution of the sfermion 
loops is enhanced, but nevertheless weighs less than 10% of the total one-loop correction. 
As the correction is dominated by fermion loops and is rather independent of for 
\/s > 500 GeV, we think that it can be approximated by a factor correction to the 
tree-level cross section. Such an approximation would be most useful for including initial 
state radiation (ISR) and beamstrahlung in an efficient way. 
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